Temperature Dependence of the Magnetic Penetration Depth and Nodal Gap 
Structure of UPts from Small Angle Neutron Scattering 
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We report measurements of the temperature dependence of the different components of the Lon- 
don magnetic penetration depth using small angle neutron scattering from the vortex lattice in a 
high quality crystal of UPta. Deconvolution of the contributions to the neutron scattering form fac- 
tor from currents along the a* and c-axes gives direct information about the nodal structure of the 
order parameter. Our observations of linear temperature dependence at low temperatures of all com- 
ponents of the penetration depth support the assignment of E2u symmetry to the superconducting 
state of UPta . 
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The heavy fermion compound UPta has attracted sub- 
stantial interest as a paradigm for unconventional su- 
perconductivityjii^ Among a host of fascinating physi- 
cal properties, perhaps most striking is the fact that the 
H — T superconducting phase diagram has three distinct 
superconducting vortex phases shovifn in Fig. 1(a), la- 
beled A, B, and C. Experiments and theory demonstrate 
that this phase diagram can only be explained by an un- 
conventional superconducting order parameter;^ a close 
parallel to superfluid '^He. However, a complete theo- 
retical description of the superconducting state of UPta 
has not been settled, and there are several candidate the- 
ories that can account for the material's unusual phys- 
ical properties. The order parameter structure consis- 
tent vifith the widest variety of experiments is an odd- 
parity, spin-triplet, /-wave orbital state of Ei2u symme- 
tryj^i^ With some success, comparisons with experiment 
have also been made for an even-parity, spin-singlet, d- 
wave orbital state of Eig symmetryi^ Both of these or- 
der parameters break time reversal symmetry in the low 
temperature B phase. A more recent proposal^ for an 
odd-parity, spin-triplet, /-wave model with Ei„ orbital 
symmetry, is time reversal symmetric in the B phase, in 
contrast to the other two models. 

All of these order parameters have nodes in the super- 
conducting energy gap, each with different nodal struc- 
ture in the three vortex phases. Consequently, it is of 
particular importance to explore physical properties that 
are directly linked to this nodal structure. In this Letter 
we report our constraints on order parameter symmetry 
based on small angle neutron scattering (SANS) mea- 
surements from the vortex lattice (VL) in the supercon- 
ducting mixed state. We have measured the temperature 
dependence of the components of the London penetration 
depth, Aj;(T), that probe the gap nodal structure along 
the principal directions of the crystalloid We compare our 
results for (T) with those from other methods including 
ac-susccptibility2iiiS and muon spin rotation (yitSR) liiii^ 
In our work we find a linear temperature dependence of A 



for currents parallel to the a* and c-axes consistent with 
an E2u order parameter. 

Our sample consists of a high-quality, 15 g single crys- 
tal (RRR > 600), cut into two pieces, and is described by 
Gannon et alJ^ The UPts crystals were co-aligned, fixed 
with silver epoxy to a copper cold finger, and mounted 
to the mixing chamber of a dilution refrigerator with the 
crystal a-axis vertical and the c and a*-axes in the hor- 
izontal scattering plane. Rotation of the dilution insert 
allowed easy reorientation of the a* or c-axes to be par- 
allel to the magnetic field and neutron beam inside a 
horizontal superconducting magnet on the SANS-I and 
SANS-II beamlines at the Paul Scherrer Institut in Vil- 
ligen, Switzerland. For measurements on SANS-I, the 
neutron wavelength was 6 A with 11 m of coUimation. 
For measurements on SANS-II, 9 A neutrons were used 
with 6 m of collimation. 

Fig. lUb) shows an example of a SANS diffraction pat- 
tern from UPt3, with magnetic field H = 0.3 T parallel 
to the a*-axis. As UPta has relatively long penetration 
depths, only the first order Bragg reflections are evident 
in this image which is a superposition from rocking angles 
about the horizontal where the Bragg condition was sat- 
isfied for diffraction spots above the beam center. For all 
data discussed in this letter, background scattering mea- 
sured in zero applied field was subtracted. The diffrac- 
tion pattern shown in Fig. [TJb) is of a distorted hexago- 
nal VL, similar to, but more anisotropic than, previous 
SANS measurements in this orientationJ^ii^ By symme- 
try, there are four additional peaks indicated by red cir- 
cles that were not imaged in the hexagonal domain in 
Fig. lUb). The symmetry of the diffraction pattern is the 
same as that of the real space VL, rotated by 90 degrees 
with a rescaling of the axes. The distortion of the VL 
from a perfect hexagon is a result of penetration depth 
anisotropy in the plane perpendicular to a*. When a 0.2 
T field is applied parallel to the c-axis, a perfect hexago- 
nal VL is seen, in agreement with previous measurements 
in that orientation at a similar fieldJ^ 



2 




(b) 



200 400 600 

r(mK) 



-0.010 0.010 

qx(A-^) 



FIG. 1. (color online) (a) A schematic of the phase diagram 
for the three vortex phases of UPts for H\\a* . (b) An example 
of a diffraction pattern, measured at « 50 mK for H = 0.3 
T with H\\a*. The Bragg condition was only satisfied for 
diffraction peaks above the beam center. By symmetry, there 
are two peaks below the beam center and 2 peaks on the 
horizontal, indicated by red circles. The UPta crystal axes 
are shown as white arrows. 
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FIG. 2. (color online) The opening angle 2a, defined in the 
inset, shown as a function of applied magnetic field along the 
a*-axis measured at ~ 50 mK. The B to C-phase transition 
is given by the vertical dashed line at iJ" = 0.6 T. Red lines 
are guides to the eye. 



Fig. [2] shows the opening angle 2a of the VL, defined 
in the inset to Fig. [2l as a function of applied magnetic 
field for H\\a*. For the data shown here the VL was 
prepared by reducing the magnetic field from above Hc2 ■ 
Then, a damped oscillation with initial magnitude 0.02 
T was performed around the desired field until the fi- 
nal measurement field was reached. The motivation for 
preparing the VL with this field history was to produce 
the superconducting state with and an equilibrium order 
parameter orientation and to ensure that the VL was in 
the ground state.— 

Our 2a data can be best described as having a linear 
field dependence in the B-phase that becomes field in- 
dependent in the C-phase, where the B to C-transition 
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FIG. 3. (color online) Rocking curves measured at ~ 50 mK 
for H = 0.3 T (black circles) and 0.85 T (purple squares) 
with fields along the a*-axis. Fits to each rocking curve are 
gaussian and the axis is shifted to the curve center, 0o. 



occurs between H = 0.5 and 0.6 T in this field orien- 
tation. Field dependence to the opening angle indicates 
that there are deviations from the London theory asso- 
ciated with non-local correctionSi '^i^^'^^ A change in be- 
havior of the opening angle at the B to C-transition was 
also reported by Yaron et alt^ 

The intensity in a diffraction peak is related to the 
Fourier transform of the local field variations from the 
real space VL. The scattered intensity in a reflection, 
/((/)) is shown in Fig. [3] as a function of rocking angle (j) 
for two different magnetic fields. The scattered intensity 
decreases as magnetic field is increased since the propor- 
tionate increase in vortex density reduces the contrast in 
local field spatial variations. We measured the first or- 
der Fourier component of the diffraction, called the form 
factor expressed as. 
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The form factor is calculated from the reflectivity, i?, 
equal to the integrated intensity, multiplied by cos a (the 
Lorentz factor), divided by the incident neutron flux. In 
Eq. 1, $o=2.07xlO^T-A2 is the magnetic flux quantum; q 
is the magnitude of the scattering vector of the reflection 
being measured; the gyromagnetic ratio of the neutron 
is 7 = 1.91; A„ is the incident neutron wavelength; and 
t is the effective sample thickness which we have taken 
to be 3.9 mm ~ the equivalent thickness of a uniform 
sample with the same width, height, and volume as our 
sample. The held dependence of our measurements of 
rocking curve widths, such as those shown in Fig. |31 do 
not show the sudden broadening at the B-C phase transi- 
tion reported by Yaron et alr^ All of our rocking curves 
are approximately 20% broader than the resolution limit 
for our experiments. We also do not see a change in slope 
of the field dependence of |/ii| at the B-C transition as re- 
ported earlier It is likely that absence of these effects 
can be attributed to the higher quality of our crystal 
and the oscillatory field procedure which we have used 
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to overcome flux pinning. 

In London theory, the form factor is related to the ma- 
terial properties through the magnetic penetration depth 
A. For an isotropic superconductor the form factor is 
given by, 



B 



1 + X^q^ 



(2) 



where ^ is the superconducting coherence length and c is 
a constant, typically taken to be 5. The fractional part of 
Eq. [2] comes directly from the London equations. -A The 
exponential factor is a correction to the London theory to 
account for the non-zero extent of the vortex coresi^ For 
an anisotropic superconductor the form factor can be ex- 
pressed in terms of the principal values of the penetration 
depth Xi corresponding to currents flowing along each of 
the principal directions of the crystal. These components 
of the penetration depth are related to corresponding di- 
agonal components of the quasi-particle mass tensor m^ 
in much the same way as for the isotropic London the- 
ory, Xf oc mi, providing a direction-specific measure of 
the low lying excitations in the superconducting state 
sensitive to gap nodes»i 

For uniaxial anisotropy, as for UPta, Ai = A2 7^ A3 
where i ~ 3 corresponds to the c-axis. The form factor 
for field along the a or a*-axis becomes;^ 
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If there is no variation in the VL geometry as a function 
of temperature - a VL phase transition or temperature 
driven disorder, for example - then the temperature de- 
pendence of the form factor given by Eq. [5] and [3] reflects 
the temperature dependence of A^. Strictly speaking, the 
coherence length ^ is also anisotropic and should be taken 
into account in the exponential factor of Eq. [31 However, 
the anisotropy in the coherence length is minimal as es- 
timated from Hc2 measurements^ and its temperature 
dependence is sufficiently weak at low temperatures, that 
the effect of anisotropy is negligible. 

We have made measurements of the temperature de- 
pendence of the VL scattering for magnetic fields along 
both the crystal c and a*-axes with the field reduced 
from above Hc2, after which the only field changes were 
damped oscillations before measurement at every tem- 
perature. Rocking curves were obtained for each orien- 
tation at base temperature and at intermediate tempera- 
tures to determine that there was no broadening as tem- 
perature was varied. The magnet and sample were ro- 
tated to the center of the rocking curve and the scattered 
intensity I{T) was measured "rocked-on" as a function 
of temperature. As an example of our raw data, we show 
I{T) at two fields parallel to the a*-axis in Fig. |4l The 
diffracted intensity does not approach a constant at the 
lowest temperatures for either field orientation, indicat- 
ing the presence of nodes in the superconducting gap. 



14 

."1 






12 - 


c 



10 - 


E 




■d 
.1—' 


8 - 


JO 


6 - 








CO 


4 - 









2 - 








- 








H\ 



I 



0.2 T 
0.4 T 

Rocking curve 
measured 



0.1 0.2 



0.3 0.4 
r(K) 



0.5 0.6 



FIG. 4. (color online) The rocked-on scattered intensity as a 
function of applied magnetic field for fields along the a* -axis 
aX H = 0.2 T (red circles) and 0.4 T (blue triangles). Solid 
red arrow indicates Tc at 0.2 T. Dashed blue arrow indicates 
Tc at 0.4 T. 
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FIG. 5. (color online) (a) The opening angle 2a and (b) the 
scattering vector g as a function of temperature. For each 
panel, data is shown for fields along the a*-axis at 0.2 T (red 
circles) and 0.4 T (blue triangles). For H\\c at 0.2 T the 
data (green squares) are within error bars of being a perfectly 
hexagonal VL. Solid lines show the average for each data set. 
Dashed lines in panel (b) show q calculated using the average 
values of 2a form panel (a), assuming B=H and single flux 
quantization. 



The opening angle and the scattering vector taken di- 
rectly from the diffraction pattern are both temperature 
independent as shown in Fig. [5l[a,b). Using their aver- 
age values, we calculated \hi\ from Eq. [T] at H = 0.2 
T for each temperature. Since the penetration depth is 
isotropic in the plane perpendicular to c, we used Eq. [2] 
and our results for \hi\ with H\\c to find Ai(T), taking 
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FIG. 6. (color online) The temperature dependence of the 
penetration depth in UPts at 0.2 T. Green squares are for 
Ai with currents along a and red circles for A3 with currents 
along c . Solid lines are fits to the data from base temperature 
up to T/Tc = 0.3. Error bars for A3 are approximately the 
size of the data points. 

c = i and ^=120 A. We then calculated A3(T) from Eq.E] 
using our values for Ai(T) and our \hi\ values for H\\a*, 
with c = I and ^=110 A. The results of this analysis are 
shown in Fig. [HI 

The limiting, low temperature region of our results 
appears to be for T/Tc < 0.3, where both Ai and A3 
are linear. In Fig. [6] we show our fits to the data in 
this region. Our extrapolation to zero temperature gives 
Ai(0) = 6,830±208 Aand A3(0) = 3,990±43A. There is 
evidence for A3 to deviate from linearity for T/Tc > 0.3, 
in contrast to Ai. 

On comparing our results with earlier work, wc find 
that interpretation of ac-susceptibility measurements re- 
quire an analysis of the real and imaginary parts of the 
electromagnetic response from which extraction of the 
penetration depth is not trivial and is necessarily sensi- 
tive to surface qualityi^iiS Signorc et al.^ reported a lin- 
ear temperature behavior which could not be associated 
with any specific component of the penetration depth. 
An early /xSR investigation by Broholm et aL— found a 
penetration depth anisotropy much too small to be con- 
sistent with other observations of the superconducting 
statei^ii^ii^ In a later ^SR study, Yaouanc et ali^ found 
Ai(0) = 6, 040 ± 130 A and A3(0) = 4, 260 ± 150 A with 
H = 0.018 T, qualitatively consistent with what wc re- 
port here. 

Our measured penetration depth anisotropy in the low 
temperature limit is A1/A3 = 1.71 ± 0.07 ai H — 0.2 
T, giving a quasiparticle mass anisotropy of mi /my, ~ 
2.94 ± 0.24 in good agreement with the anisotropy of 
2.93 ± 0.3 that wc calculate from the VL opening angle 
ai H = 0.2 T using the relation tana = (7713/37711)^/^ 
taken from the London theory.— In earlier SANS mea- 
surements, Kleiman et alr^ and Yaron et alr^ report 
analysis of the opening angle and mass anisotropy. Us- 
ing their published values of the opening angle at _ff = 0.2 



T, we infer their mass anisotropics to be 2.66 and 2.15, 
respectively. It is difficult to make a reliable compari- 
son of our results for A1.3 with these two experiments, 
as they assume field independent penetration depths, in- 
consistent with measurements of the field dependence of 
the opening angle, and which might not be theoretically 
justified.— 

To interpret our data in terms of the pairing symme- 
try of UPt3, we provide a brief discussion of the nodal 
structures of the superconducting gap for the various can- 
didate theories. In the low field and low temperature 
B-phase where all of the data shown in Fig. |6] were mea- 
sured, the three predominate pairing models discussed 
earlier have three different nodal structures. For the E2„ 
model)^ there are point nodes at the poles of the Fermi 
surface which open with quadratic wave-vector disper- 
sion, and a line node around the equator of the Fermi sur- 
face that opens with linear dispersion. The Ei^ model^ 
also has point nodes at the poles, however these nodes 
open linearly Similar to E2„, the Ei^ model also has 
a line node around the equator that opens linearly. The 
Ei„ model^ has a somewhat more complicated gap struc- 
ture in the B-phase, with point nodes at the poles that 
have a non-trivial dispersion and an antinode at the equa- 
tor. 

Our measurements of Ai are sensitive to the nodal 
structure at the equator, while those of A3 are sensitive to 
the nodal structure at the poles. For a point node with 
quadratic dispersion, a linear temperature dependence 
of A is expected in the low temperature limit, while for 
a point node with linear dispersion, a T^ temperature 
dependence of A is expected. Our observed linear depen- 
dence of A3 is evidence for quadratic point nodes. For 
a line node with linear dispersion, a linear temperature 
dependence in A is predicted. And for an antinode along 
the equator, power-law temperature dependence should 
not occur Our observed linear dependence of Ai is ev- 
idence for a line node at the equator. Our conclusions of 
quadratic point nodes and an equatorial line node in the 
B-phase gap support E2M pairing symmetry. 

In summary, we have observed linear temperature 
dependence for both Ai and A3 at low temperatures, 
providing strong evidence for an equatorial line node 
and quadratic polar point nodes in the energy gap of 
UPt3 . Our results support the assignment of E2U pairing 
symmetry and are consistent with the conclusions from 
some other recent experiments j 1 ^"^ However, there 
arc still a number of unresolved questions concerning 
the nature of the order parameter in UPt3 including 
the importance of spin-orbit interaction and the spin 
structure of the order parameter,— i^ii^i^^i^ as well as a 
direct observation of broken time reversal symmetry in 
the superconducting B-phase. 
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